Weak ferromagnetism and magnetoelectric coupling in bismuth ferrite 
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We analyze the coupling between the ferroelectric and magnetic order parameters in the magneto- 
electric multiferroic BiFeOs using density functional theory within the local spin density approxima- 
tion and the LSDA+U method. We show that weak ferromagnetism of the Dzyaloshinskii-Moriya 
type occurs in this material, and we analyze the coupling between the resulting magnetization and 
the structural distortions. We explore the possibility of electric-field-induced magnetization reversal 
and show that, although it is unlikely to be realized in BiFeOa, it is not in general impossible. 
Finally we outline the conditions that must be fulfilled to achieve switching of the magnetization 
using an electric field. 
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There has been increasing recent interest in magneto- 
electric multiferroics 0, 0, 0, 0, 1 which are materials 
that show spontaneous magnetic and electric ordering in 
the same phase. In addition to the fascinating physics 
resulting from the independent existence of two or more 
ferroic order parameters in one material , the coupling 
between magnetic and electric degrees of freedom gives 
rise to additional phenomena. The linear and quadratic 
magnetoelectric (ME) effects, in which a magnetization 
linear or quadratic in the applied field strength is in- 
duced by an electric field (or an electric polarization is 
induced by a magnetic field), are already well established 
. Recently, more complex coupling scenarios have been 
investigated. Examples are the coupling of the antiferro- 
magnetic and ferroelectric domains in hexagonal YM11O3 
or the large magnetocapacitance near the ferromag- 
netic Curie temperature in ferroelectric BiMn03 . Es- 
pecially interesting are scenarios where the direction of 
the magnetization or electric polarization can be modi- 
fied by an electric or magnetic field respectively. Such 
a coupling would open up entirely new possibilities in 
data storage technologies, such as ferroelectric memory 
elements that could be read out nondestructively via the 
accompanying magnetization. Some progress has been 
made in this direction. Recently, the small (0.08/iC/cm 2 ) 
electric polarization in perovskite TbMnC>3 was rotated 
by 90° using a magnetic field at low temperatures (~10- 
20K) 0. Conversely, early work on nickel-iodine boracite 
showed that, below ^60K, reversal of the spontaneous 
electric polarization rotates the magnetization by 90° in- 
dicating that the axis of the magnetization, but not its 
sense, can be controlled by an electric field. In fact, it 
was believed 8] that electric-field-induced 180° switching 
of the magnetization should be impossible, because a re- 
versal of the magnetization corresponds to the operation 
of time-inversion whereas the electric field is invariant 
under this operation. In this work we show that such be- 
havior is not generally impossible by using multiferroic 
bismuth ferrite, BiFeOs, as a test case to analyze the 



coupling between magnetism and ferroelectricity. 
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FIG. 1: Schematic view of the R3c structure built up from 
two cubic perovskite unit cells. The cations are displaced 
along the [111] direction relative to the anions and the oxygen 
octahedra rotate with alternating sense around the [111] axis. 

BiFeOs has long been known in its bulk form to be 
an antiferromagnetic, ferroelectric multiferroic 0, llfj| . 
with antiferromagnetic Neel temperature Tn ^643K, and 
ferroelectric Curie temperature Tq ~T103K. It has a 
rhombohedrally distorted perovskite structure with space 
group i?3c [ll], ( see Fig. QJ. The Fe magnetic mo- 
ments are coupled ferromagnetically within the pseudo- 
cubic (111) planes and antiferromagnetically between 
adjacent planes (so-called G-type antiferromagnetic or- 
der). If the magnetic moments are oriented perpendic- 
ular to the [111] direction, the symmetry also permits a 
canting of the antiferromagnetic sublattices resulting in 
a macroscopic magnetization; so called weak ferromag- 
netism [Hi . Il4| . However, superimposed on the anti- 
ferromagnetic ordering, there is a spiral spin structure 
in which the antiferromagnetic axis rotates through the 
crystal with an incommensurate long-wavelength period 
of -620A 0. This spiral spin structure leads to a can- 
cellation of the macroscopic magnetization and also in- 
hibits the observation of the linear ME effect ^(|. How- 
ever, a significant magnetization (^1/xb per unit cell), 
as well as a strong ME coupling, have been reported re- 
cently in high quality epitaxial thin films 0. This ob- 
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servation suggests that the spiral spin structure is sup- 
pressed in thin films, perhaps due to epitaxial constraints 
or enhanced anisotropy. Since these epitaxial films also 
show large electric polarization (^50-60^C/cm 2 ) they 
are promising candidate materials for ME device appli- 
cations. In this work we present results of first principles 
calculations of the magnetic properties of BiFeC>3 and an- 
alyze the coupling between the magnetic and ferroelectric 
properties in this material. 

Our approach is based on density functional theory 
(DFT), see e.g. an d we use two different im- 

plementations to cross-check our results; the projector- 
augmented plane-wave (PAW) method |l8( implemented 
in the Vienna Ab-initio Simulation Package (VASP) 0, 
Eflj and the linear muffin-tin orbital method in the atomic 
sphere approximation (LMTO-ASA) j^ij, extended for 
the treatment of non-collinear spin-configurations and 
spin-orbit coupling j^,!^. These two features, which 
are also implemented in the VASP code, are often omit- 
ted in DFT calculations, but are essential for our in- 
vestigations. Except where explicitly stated, we use the 
crystal structure obtained by optimizing the atomic po- 
sitions within the experimentally observed R3c symme- 
try (see Fig. 0. Our calculated structural parameters 
are identical (within the usual numerical accuracy) with 
those given in Q and agree well with the experimentally 
observed structure To reproduce the situation in 

the thin films we suppress the spiral spin structure in 
our calculations. We also use two different treatments 
of the exchange-correlation functional; the standard lo- 
cal spin-density approximation (LSDA, see an d the 
LSDA+U method The LSDA+U method introduces 
two parameters into the treatment of the Fe d states, the 
Hubbard parameter, U, and the exchange interaction, J. 
We use a value of J = leV and treat U as a free pa- 
rameter, varying it from leV to 7eV, while keeping the 
structure fixed to that calculated within the LSDA. For 
the PAW calculations we use potentials with 15 valence 
electrons for Bi (5d 10 6s 2 6p 3 ), 14 for Fe (3p 6 3d 6 4s 2 ), and 
6 for O (2s 2 2p 4 ) and a plane wave cutoff of 400eV. For the 
Brillouin-zone inte grat ions we use a 5x5x5 Monkhorst- 
Pack k-point mesh l25l and the tetrahedron method with 
Blochl correction [U (both PAW and LMTO). These 
values result in good convergence for all quantities under 
consideration. 

Since the nature of weak ferromagnetism in BiFe03 
is not well established due to the presence of the spi- 
ral spin structure in the bulk, we first investigate the 
occurrence and origin of weak ferromagnetism, which is 
intimately connected with the symmetry of the system 
[l3|. In BiFe03 it can only occur if the sublattice mag- 
netizations are oriented in the (III) plane so that the 
symmetry is reduced to the magnetic space group Bb or 
Bb' which (apart from the primitive translations) con- 
tains only one glide plane. In this case a canting of the 
magnetic sublattices does not lead to a further reduc- 



tion in symmetry and weak ferromagnetism can occur. 
Therefore we first determine the preferred orientation 
of the sublattice magnetizations. We do this by calcu- 
lating the energy difference between the arrangements 
with (i) magnetic moments aligned parallel/antiparallel 
to the [111] direction and (ii) magnetic moments oriented 
within the (111) plane. Both our methods result in an 
LSDA energy difference of about 2meV, with a preferred 
orientation of the magnetic moments within the (111) 
plane. This arrangement is compatible with the exis- 
tence of weak ferromagnetism. Within the (111) plane, 
orientations of the sublattice magnetizations parallel or 
perpendicular to the glide plane are energetically equiv- 
alent. The anisotropy energy is reduced to more realistic 
values within the LSDA+U method, but the easy mag- 
netization orientation is unchanged. 

The anisotropy calculations show only that weak ferro- 
magnetism is symmetry-allowed, not that it will actually 
occur. Therefore we next calculate the magnitude of the 
effect. We initiate our calculations to a homogeneous 
and collinear spin arrangement with the magnetic mo- 
ments oriented in the (111) plane (along either the x or 
y-axis in our coordinate system, see Fig. , then let the 
magnetic moments relax freely within the self-consistency 
cycle. The magnetic moments then cant away from the 
collinear direction (while remaining in the (111) plane) 
by an angle of about 1° (LSDA, Fig.Eb). This leads to 
a small but measurable magnetization of approximately 
O.I^b per unit cell. LSDA+U calculations give the same 
qualitative results but with slightly smaller magnetiza- 
tions. This value is smaller than that reported in Q but 
agrees well with more recent measurements |27j . 



(a) 





FIG. 2: (a): Starting configuration of our calculation. The 
magnetic moments, Mf c i and Mf C 2, of the two iron atoms in 
the unit cell are oriented antiferromagnetically and collinearly 
in the (111) plane, allowing weak ferromagnetism by symme- 
try, (b): Calculated magnetic structure including the spin- 
orbit interaction: The two iron magnetic moments rotate in 
the (111) plane so that there is a resulting spontaneous mag- 
netization, M. 

According to Dzyaloshinskii and Moriya (DM) [l3L Hif 
the canting of the magnetic sublattices is caused by an 
antisymmetric spin coupling, the so-called DM interac- 
tion, which is due to the combined action of exchange 
interaction and spin-orbit coupling. Indeed, if we neglect 
the spin-orbit interaction in our calculations the mag- 
netic moments remain collinear and there is no macro- 
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scopic magnetization. The DM interaction has the form 

E DM = -^D • (Mfci x M Fc2 ) = -D • (L x M) , (1) 

where D is a coupling vector analogous to the Heisenberg 
exchange constant J in the usual symmetric exchange 
interaction. The antiferromagnetic vector L = Mpoi — 
Mp 2 is defined as the difference of the two sublattice 
magnetizations, and M = Mp c i + Mp G 2 is the resulting 
magnetization. From the form of -Edm it is clear that, for 
constant D and fixed orientation of L, the canting of the 
magnetic sublattices always occurs such that D (required 
by symmetry to be oriented along the [111] axis), L, and 
M build up a right-handed system. Indeed, if we start 
our calculation with a spin configuration in which the 
magnetic moments are canted in the wrong direction (so 
that D, L, and M make up a left-handed system) the 
moments relax back into the right-handed configuration 
during the iteration process. Therefore, for a particular 
orientation of D and L, only one canting direction lowers 
the energy relative to the collinear state. 

Next we analyze the relationship between the weak fer- 
romagnetism and the structural distortions in BiFeC>3. 
As already mentioned, weak ferromagnetism depends 
crucially on the symmetry of the system, which in turn is 
determined by the structural distortions. The ferroelec- 
tric i?3c structure of BiFeC>3 is reached from the ideal 
cubic perovskite structure by freezing in two unstable 
normal modes: (i) the polar displacements of all the an- 
ion and cation sublattices relative to each other, which 
lead to the spontaneous electric polarization, and (ii) 
an antiferrodistortive rotation of the oxygen octahcdra 
around the [111] direction with alternating sense of rota- 
tion along the [111] axis (see Fig. [I}. In terms of symme- 
try groups, the polar displacements alone would reduce 
the symmetry of the ideal perovskite structure (Pm3m) 
to the rhombohedral space group i?3m, whereas the ro- 
tation of the oxygen octahedra alone would lead to space 
group i?3c. The incorporation of both distortions gives 
the actual space group of BiFeC>3, R3c. Weak ferromag- 
netism is only allowed by symmetry in the space groups 
i?3c and i?3c, suggesting that it is related to the oxygen 
rotations rather than to the polar displacements along 
[111]. We have verified this by performing calculations 
for structures containing only one of the two distortions 
(while keeping the lattice vectors fixed to those of the 
i?3c structure). These calculations confirm that the po- 
lar R3m structure does not show weak ferromagnetism, 
whereas the non-ferroelectric i?3c does. 

To fully understand the coupling between the struc- 
tural distortions and the magnetization, we next invert 
the sense of rotation of the oxygen octahedra while keep- 
ing the polar distortions fixed. Again we start from the 
magnetic configuration shown in Fig. i.e. with the 
same orientation of L as in the previous calculations. 
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TABLE I: All possible orientation states with paral- 
lel/antiparallel direction of M. The first row gives the lost 
symmetry element that maps state 1 onto state i. 1 stands 
for space-inversion, a prime indicates time-inversion, m is the 
mirror plane parallel to the [111] direction, and 2 is a twofold 
axis perpendicular to m. The states 5-8 are the antiphase 
domains not included in Aizu's scheme (see text). The direc- 
tions of P, D, M, and L projected on the (111) plane are 
indicated by arrows. (©) indicates orientation along the 
positive (negative) [111] axis. 



In this case the magnetization direction is reversed from 
that of the original structure. Conversely, if we invert 
the polar distortion while keeping the rotational sense of 
the oxygen octahedra fixed, the magnetization direction 
is unchanged. This clearly shows that the direction of 
the DM vector D is determined by the sense of rotation 
of the oxygen octahedra surrounding the corresponding 
magnetic ions, rather than by the direction of the polar- 
ization as suggested in . 

Next we investigate the switching possibilities of the 
system from one stable orientation state to another by 
the application of an electric field. A general scheme 
for the derivation of all possible orientation states in fer- 
roic materials was developed by Aizu |29j |. who showed 
that the stable states can be constructed by applying all 
symmetry elements that are lost during the ferroic phase 
transition to an arbitrary orientation state of the ferroic 
phase. A complication arises in the case of BiFeC>3 be- 
cause the symmetry change between the nonferroic pro- 
totype phase (nonmagnetic PmZm) and the final multi- 
ferroic phase (i?3c) is not purely ferroic in nature. The 
unit cell doubling, caused by the oxygen rotations, trans- 
forms the mirror plane parallel to the [111] direction into 
a glide plane, and leads to the formation of antiphase do- 
mains which are not described by Aizu's scheme. To in- 
clude the antiphase domains one has to consider the full 
space group symmetry instead of only the point group 
symmetry; for BiFeOs this analysis leads to a total of 96 
degenerate orientation states. Here we limit our discus- 
sion to the 8 different orientation states with magnetiza- 
tion parallel/antiparallel to a fixed axis. These represent 
180° switching of the magnetization, and could be iso- 
lated in practice by lifting the in-plane degeneracy of the 
easy magnetization axis using epitaxial strain. The re- 
sulting orientation states are listed in Table [I] indicated 
by the different orientations of M, L, D and P, where P 
is the electric polarization and the "sense" of the oxygen 
rotations is indicated by D. 
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The application of a polarization-reversing electric field 
could in principle drive the system from initial state 1 
into any of the degenerate states 2, 4, 6, or 8. Of course 
in reality the system will prefer to change into the state 
separated from the initial state by the lowest energy bar- 
rier. A reversal of L is unlikely since it involves a ~180° 
rotation of the sublattice magnetizations, which is hin- 
dered by the magnetic anisotropy; in contrast the reversal 
of M requires only a small reorientation of the magnetic 
moments. This reduces the probable outcomes to either 
state 4 (in which both D and M reverse) or state 8 (in 
which D and M are unchanged). Of these, state 8, in 
which the magnetic ordering is unchanged compared to 
state 1, is the most likely, since the reversal of the oxy- 
gen rotations is energetically costly, and is not required 
by the reversal of P. It is clear, however, that the ear- 
lier argument || that the mutual invariance under timc- 
and space-inversion of electric field and magnetization in- 
hibits the possibility of electric-field induced magnetiza- 
tion reversal, does not hold. While it is certainly correct 
that it is not possible to drive the system from state 1 
to its time-conjugate state 1' using an electric field, the 
state corresponding to V is also present in all ferroelec- 
tric ferromagnets. If the energy barrier for the transition 
to 1' or to any other orientation state with reversed M 
and P is lower than the energy barrier to all other de- 
generate orientation states, then the magnetization will 
be invertable by an electric field. 

From the above discussion we can extract three con- 
ditions that must be fulfilled to achieve electric-field- 
induccd magnetization reversal in such a rhombohedrally 
distorted multiferroic perovskite: (i) the rotational and 
polar distortions must be coupled, (ii) the degeneracy be- 
tween states 1, 2, 3, 4 and 5, 6, 7, 8 must be lifted, i.e. 
parallel and antiparallel orientations of D and P must 
be inequivalent, and (iii) there should be only one easy 
magnetization axis in the (111) plane. The latter condi- 
tion can easily be achieved by straining the material in an 
appropriate way. Fulfillment of the other two conditions 
will take further exploration but we are not aware of any 
general restriction that would make this impossible. 

In summary, we have shown that BiFe03 exhibits weak 
ferromagnetism of the DM type if the spiral spin struc- 
ture is suppressed. We have also shown that the DM vec- 
tor is determined by the rotations of the oxygen octahe- 
dra rather than by the ferroelectric polarization. Finally 
we have discussed the possible magnetoelectric switch- 
ing scenarios in BiFe03, and formulated conditions that 
must be met to realize electric-field-induced magnetiza- 
tion reversal. 
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